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Abstract 

The long wavelength physics in a de Sitter region depends on the initial quantum state. While 
such long wavelength physics is under control for massive fields near the Hartle-Hawking vacuum 
state, such initial states make unnatural assumptions about initial data outside the region of causal 
contact of a local observer. We argue that a reasonable approximation to a maximum entropy 
state, one that makes minimal assumptions outside an observer's horizon volume, is one where a 
cutoff is placed on a surface bounded by timelike geodesies, just outside the horizon. For sufficiently 
early times, such a cutoff induces secular logarithmic divergences with the expansion of the region. 
For massive fields, these effects sum to finite corrections at sufficiently late times. The difference 
between the cutoff correlators and Hartle-Hawking correlators provides a measure of the theoretical 
uncertainty due to lack of knowledge of the initial state in causally disconnected regions. These 
differences are negligible for primordial inflation, but can become significant during epochs with 
very long-lived de Sitter regions, such as we may be entering now. 
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I. INTRODUCTION 



At the classical level, de Sitter spacetime appears to be stable. For conformally-coupled 
matter rigorous nonlinear stability theorems have been established that show for an open 
set of initial data the solution evolves to an asymptotically de Sitter solution at late times 

M. 

Quantum mechanically, the situation is less clear. In a quantum theory correlators do 
depend on data outside the past light cone of points. It is therefore important to study the 
effect of choice of initial state on predictions for cosmological correlators, that then in turn 
determine the spectrum of density perturbations at late times. 

There has been much confusion in the literature on this point. Broadly speaking there 
are two main camps, those who advocate using a dS invariant quantum state, to compute 
dS invariant correlators; and those who impose more physically motivated initial conditions, 
breaking the de Sitter symmetries. In the first scenario, the slow roll parameters provide the 
only means of explicit breaking of the de Sitter isometries. It has been convincingly argued 
that at least in massive scalar theories, these correlators are infrared finite |50j. It remains 
an open question whether the same is true when the quantum fluctuations of gravitons are 
included. 

In the other camp are those who have typically had in mind direct applications to inflation. 
In that case, an initial state is chosen at some finite time, breaking the de Sitter invariance. 
Needless to say, the correlators break de Sitter invariance, and in massless scalar theories, 
exhibit infrared divergences. Infrared cutoffs lead to terms in the correlators that grow 
as log a(t) to some power [31422] (here a(t) is the scale factor in the Friedmann-Robertson- 
Walker metric). In this setting there have also been computations involving gravitons at one- 
loop, which also find the loga(t) growth |23| . This has led to the suggestion that quantum 
effects may lead to a relaxation of the cosmological constant [5j . See [24J for a review article 
summarizing many of these approaches. Some also [25, 26J for some observations closely 
related to the approach discussed in the present paper. 

In this paper we will argue the main results of these two camps are mutually compat- 
ible, and that the difference between the two approaches may be viewed as a theoretical 
uncertainty in the relevant correlator. Since it is likely we do not have causal access to the 
spacetime region that would allow us to precisely determine the initial state, the best we can 
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do is perform some version of a Bayesian analysis and quantify the theoretical uncertainty 
in the initial state, and hence the derived correlators. 

If we work within the general framework of inflation, the initial state with the minimal 
number of assumptions that gives rise to our present observable universe, is a state that 
started as an approximately homogeneous region an inverse 10 14 GeV in size (see |27] for 
a more detailed discussion of the initial conditions for inflation). To explain the horizon 
problem, this region must inflate to a size of about 1 m, assuming inflation ends, and 
reheating produces a thermal gas with temperature around 10 14 GeV (this gives the famous 
60 e- foldings of inflation needed for viability). Subsequently the universe evolves according to 
the Standard Model of cosmology. Thus, at least classically, we can conclude that any state 
that respects these conditions will give rise to what we see today. According to the Bayesian 
approach, we should pick the most likely such state (or maximum entropy state). The 
most obvious such state would involve the homogeneous patch, but would be surrounded 
by a thermal gas at the Planck temperature. Of course such a state would involve large 
gravitational back-reaction, and it would be impossible to compute with it using known 
methods |5T]. In the following we will adopt a compromise that leaves one with computable 
correlators, but removes most assumptions about the initial region outside the homogeneous 
region. Therefore we propose to use an initial state with a hard infrared cutoff just outside 
the horizon at the start of inflation when the homogeneous patch begins to expand. 

After this time, the homogeneous patch expands. If the infrared cutoff was kept at a 
fixed proper length, the patch would soon exceed the size of the cutoff, and the resulting 
correlators would miss much of the physics of interest to us today. Moreover, to keep the 
cutoff at the same length, the walls of the box would have to contract faster than the speed 
of light, which signals an unphysical choice of regulator. Nevertheless, such a cutoff is needed 
for massless theories in pure dS spacetime to retain the dS invariance of the correlators. 

The more physical choice we advocate in the present work is to view the infrared cutoff 
as a kind of bubble wall. The trajectory of the wall will depend on the details of the 
bubble, but for any physical choice, at late times, the wall will asymptote to a family of 
timelike geodesies. Therefore the simplest cutoff that captures these qualitative features 
is a comoving cutoff. Such a cutoff has been studied before in the literature in 1 17. 28], 
where it is referred to as a minimal box. Thus we advocate using a comoving infrared cutoff 
placed just outside the horizon at the beginning of inflation. This explicitly breaks de Sitter 
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invariance. As we will see, this can lead to important effects at late times in correlators of 
massless fields. 

The dS invariant computations, on the other hand, represent the maximal number of as- 
sumptions about the initial state. Choosing the Euclidean vacuum across the entire inflating 
patch amounts to solving the horizon problem by hand by imposing gaussian fluctuations of 
fixed size across the entire initial slice. The chief advantage of doing this is that, at least for 
massive scalar theories (where no additional infrared cutoff is needed), the correlators are 
dS invariant, and can be much more simply treated to extract predictions for observations 
today. Moreover, as we shall see, the answers agree with the leading terms of the correlators 
with a comoving cutoff, provided inflation does not last too long. The difference between 
the two approaches should be interpreted as a measure of the theoretical uncertainty in the 
correlator due to our ignorance of the initial state at the start of inflation. 

As we shall see, the loga(t) terms in the examples we study are subleading versus the dS 
invariant results provided inflation does not last extraordinarily long. Thus predictions for 
primordial inflation are largely unchanged. 

On the other hand, we may be entering a regime of dS dominance in our present epoch. 
Assuming the graviton is the only field that experiences these loga(t) divergences, we can 
expect this quantum instability of de Sitter space to become relevant only after a very large 
number of e-foldings. We note this instability is predicted from a unitary model of quantum 
gravity based on embedding dS regions in asymptotically AdS spacetimes, dual to conformal 
field theories PEED]. 

II. SCALAR FIELD THEORY USING THE IN-IN FORMULATION 

To illustrate the interpretation of amplitudes mentioned above, we take examples from 
4 scalar field theory. Our goal is to exhibit the one-loop corrections computed using the 
comoving infrared cutoff (our approximate maximum entropy initial state) and compare 
them to computations done without an infrared cutoff, choosing the usual Euclidean vacuum 
state. It is also necessary to specify an ultraviolet cutoff to regulate the one-loop integrals. 
We choose an ultraviolet cutoff motivated by local effective field theory, simply a cutoff at 
fixed proper momentum in both cases. While these or closely related results have already 
appeared in the literature [16[ EH EI], we review the derivation in some detail to establish 
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a consistent notation and collect all the relevant results together. 
We work with the Lagrangian 



C 



-fj 



^d^dy - ^m 2 (f) 2 - I *C with 



5C 



where 5C contains the counter terms. The de Sitter metric is written as 



ds 2 = -dt 2 + a(t) 2 dx 2 = a(T) 2 (-dT 2 + dx 2 ), 



(1) 



with a(t) = e Ht , and conformal time is defined as r = —H 1 a(t) 1 . We use the notation 
x = (x, t). 



A. Mode expansions 



To set up the perturbative expansion, we begin with the expansion of the free field in 
modes labeled by comoving wavevectors k 

d 3 k 



J< >(af,r) 



<Pk{r) 



(2tt) 3 

y/—TTT 



e ifc %(r)^ + e-*^(r)4 



with 



2a t) 



(2) 
(3) 



Here is the Hankel function of the first kind and v 2 = 9/4 — m 2 /H 2 . The satisfy 
standard canonical commutation relations. The Bunch-Davies vacuum is defined by o;r|0) = 
0, for all k. 



B. In-In formalism 



We will work with the in-in formalism, where the goal is to compute expectation val- 
ues of time-ordered products of Heisenberg field operators with respect to our fixed initial 
state. This differs from the usual path integral approach, which computes transition ampli- 
tudes. A nice review of the in-in, or Schwinger-Keldesh approach can be found in |32j . The 
perturbative expansion can be derived using interaction picture methods as usual 

(0|T ((KxuhMx^h)) |0> = (0|T Co (e- i ^^(*') d *>W(f 1 ,t 1 )0 (o) (^,t 2 )) |0> 
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Figure 1: This shows the Keldesh time-ordering contour Co = C+ U C_, for the case t 2 > t\. The 
initial state is specified at time to- 

with the new feature being the appearance of the Schwinger-Keldesh time-ordering operator 
Tc which accomplishes the task of time evolving the amplitude from the initial time to, 
to max(ti,t 2 ) and then back again to to so that the expectation value is computed. The 
contour Co is shown in figure [Tj The operator insertions at t\ and t 2 appear on the upper 
contour C + . Here H{ nt is the interacting part of the Hamiltonian obtained from 

When the exponential is expanded in powers of A we can apply Wick's theorem to break 
the expectation value into integrals of products of the following elementary time-ordered 
two-point functions: 

G T (xi,x 2 ) = -iCr[0(xi,ti)0(x a ,£ 2 )]), h,t 2 E C + , 
G < (x 1 ,x 2 ) = i(0(x 2 ,t 2 )0(xi,ti)), t\ g C + , t 2 eC-, 
G > (x 1 ,x 2 ) = -i((j)(x 1 ,t 1 )(j)(x2,t2)), h e C_, t 2 eC + , 
k G T (x 1 ,x 2 ) = -i{fy)(xi,t 1 )^(x 2 ,t 2 )]), h,t 2 e C_. 
We then have the relations 

G T (x 1 ,x 2 ) = 0(ti - t 2 )G < (xi,x 2 ) + Q(t 2 - t 1 )G > (xi,x 2 ), 
G T (x u x 2 ) = - t 2 )G > (x l ,x 2 ) + 0(t 2 - t 1 )G < (x 1 , x 2 ). 

The contour Co can be written as a single K by writing 

\ M^,t), t e C+, 

<j){x,t) = < 

I <j)_{x,t), t e CL, 

and the two-point functions can be combined into a 2x2 matrix G as 

4(T[(j) + {x u ti)(j) + {x 2 ,t 2 )\) z(0 + (xi,ti)0_(f 2 ,*2)) 

i, ti)0+(f 2) ^)) *i)0_(f 2) t 2 )]) 

It is convenient to work in a different basis in field space by defining 
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1 X ~~ -ia(Ti) 4 5(ri-T 2 )5, 7 
Figure 2: The Feynman rules. 



and in this basis the matrix of correlators becomes 



G K = RGR T = | | , with 



iF G R 
G A 



iF(xi,x 2 ) = ^ (G > (x 1 ,x 2 ) + G < (x 1 ,x 2 )) , 
G R ( Xl ,x 2 ) = Q(h - t 2 ) (G<( Xl , x 2 ) - G>( Xl , x 2 )) , 
G A (x u x 2 ) = e(t 2 - h) (G > (x 1 ,x 2 ) - G < (x 1 ,x 2 )) . 

C. The Feynman rules and the one-loop diagrams 

In the (0c, 0A)-basis the Lagrangian becomes 

C = £(0 + ) - C{4>_) = y/^g (-d^ c d^ A - m 2 c <M - ^(<p c <Pl + 40 3 c 0a) 

5x (J, j3 1/1^3 



and we read off the Feynman rules of figure [2] 

The Feynman rules related to the counter-terms 5z and 8\ yield additional powers of the 
scale factor a(t) and coupling A respectively, and will not contribute to a late time, one-loop 
computation. Hence we only include the Feynman rule corresponding to 5 m . Since the equal 
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Figure 3: The diagrams contributing to F(xi,x 2 ) at order A. 

time propagator G R (x, x) vanishes by construction, any non-vanishing one-loop diagram has 
to have a solid line inside the loop. Thus by the above Feynman rules the only contributions 
at one-loop level are given in figure |3j According to the Feynman rules, these diagrams 
contribute 

/d 3 k 
e ^-5 2 ) Fii(A . )Ti)T2)) with (4) 

F 1L (k,T 1 ,r 2 ) = J dr v F(k J r u r v )(-tG R (k,T 2 ,T v ))(L(r v ) + C(r v )), (5) 
L(t„) = -i\a(r v ) 4 ——-F(p,T v ,T v ), C(r v ) = -ia(T v ) 4 S r , 



where L(r v ) is the contribution from the amputated one- loop diagram, and C(r„) comes 
from the counter-term diagram. 

D. The loop contribution 

As explained in the introduction, we consider a comoving infrared cut-off Air and a 
physical ultraviolet cutoff Auv- The contribution from the loop integral can then be written 

as 



"A uv a(r v ) 

L(r„) = -i\a(T v f I / +/ I 7 ^2p 2 F(p,T v ,r v ), (6) 



IR 



'''' > 

where we split the integral into parts corresponding to modes inside and outside the horizon 
at time t v . In these regions the propagators can be expanded using the mode functions (J3| 

as 

[ F(ife,Ti,r 2 ) « £s(k 2 TiT 2 ) e , 
Outside: kn < 1 < 2k ' (7) 

G R (k, Tll r 2 ) « 0(n - r 2 )f [rf-VI - rf^] , 

F(A;,r 1 ,r 2 )^^cosA;(r 1 -r 2 ), 
Inside: ArT* > 1 < (8) 

G^n.^weCn-^^pisiii^ri-ra), 
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where e = m 2 /(3H 2 ) is the mass parameter. Inserting the expansions into the loop integral 
@ we get 

"A uv a(T v ) 



u 



-iX 



_2c 



(2n)m^ y JAlR 



dp p- 1+2e + T 2 V 



dp p 



(27r) 2 i7 2 r4 



2e 



1 Muv 

2 \ H 



(9) 



Note that the integrals are dominated by the IR and UV regions as opposed to the horizon 
\pr v \ ~ 1, and the approximations (7l8) can be used. In order to cancel the ultraviolet 
divergence we fix the counter-term coefficient S m to be 



XH 2 



2(2tt) 



A 



uv 



H 



— y 

h) 



where /i is the renormalization scale. This leads to the contribution 



^ + - ^ + 1 (f) 2 ) - ^v (T ,. do) 

E. The propagator at late times 

We then wish to analyze the late time dependence of the propagator Q. Inserting the 



loop contribution (10) we get 

F iL (xi,x 2 ) = -^-^ J e *(*"*) J -^F{k^r v )G R {k^r v )V{r v ). (11) 

The integration region naturally splits into components depending on whether mode k is 
inside or outside the horizon, as indicated in figure |4j In the following we will take T\ < t 2 for 
convenience, and also assume that they are in the same cosmological period, t\ ~ r 2 ~ T now . 
The limit r 2 3> T\ is also of interest, and we will treat it separately in the following. 

The IR modes: Region I consists of modes outside the horizon. Using the expansions 
(J7J) we find the momentum space propagator to be 

Fi L (k,T 1 ,T 2 ) = -—J ——^-[T 2 r v -r v r 2 ]V(r v ), 
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Figure 4: Plot of the integration region in (r, /c)-space. 



where the lower limit of integration is given by Tx = Max(ro, —k 1 ), see figure [4j This can 
be evaluated analytically, yielding 



1 fl i 



XH 2 (k 2 nT 2 y (A ir t 2 



1 2c 



(2tt) 2 6k 3 



2e 
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2e 
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2e 



3-2e 



V 



3 - 2e 



- log 
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T2 



XH 2 (fc 2 nr 2 ) 
(2n) 2 6/c 3 



-£>(rx). 



:i2) 



Transforming back to position space we have 
d 3 k 



F( L (xx,x 2 ) 



(27T) 



e iH^2) F i L{kjTuT2] 



(27T) 



sin(/cAa;) 



k 2 Fl L (k,n,T 2 ) 



dkk 2 Fl L (k,T X ,T 2 ){l + 0((jfeAx) 



where Ax = |xi — 22 1, and in the last step we took the physical separation between X\ and 



X2 to be much less than the Hubble distance. Inserting (12) we find 



(2vr) 4 

XH 2 {t\T2) 

3(2tt) 4 



,2c 



dk— D(max(r , —k 1 )) 
3k 



D(to) 



dkk' 



-l+2e 



+ 



dkk- 1+2t D(-k- u 



(13) 



The first integral is over modes outside the horizon already at r , and vanishes if one takes 
the cut-off to coincide with the initial horizon. The integrals can be computed analytically 
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to yield for late times 

Fl L (xi,x 2 ) = 



XH 2 



12(27r) 4 e 



T\T 2 



lot 



To 
T 2 



+ 



(1 - (A IR r ) 2e ) 



+ (^)! log ^ + l fn 

e t 2 2e \t 2 



h) 



+ 



1 



1 




(14) 



The first two terms vanish for late times, while the last one asymptotes to a constant value 
that is independent of the infrared cutoff Air. 



The last term in (14) looks worrisome, as it is divergent in the limit where one takes the 



leg r 2 to future infinity while keeping the other leg T\ fixed. This is due to our assumption 
that the legs are in the same cosmological period, which no longer holds when r 2 is taken to 
the far-future. The divergence arises because the horizons at t 2 and n no longer coincide 
in this limit, but rather modes exit the T\ horizon earlier than the r 2 horizon. In the limit 
r 2 — > with 7~i fixed, with the correct expansions (7|8), the propagator (11) becomes 



F\ L {x x ,x 2 ) 



XH 2 (nr 2 y 



(2e) 2 2 



-V 

hj 



+ 



log 



n 

r 2 



(cos 1 + 0(e)) 



(2tt) 2 3 
which is well-behaved for t 2 — > 0. 

The intermediate region: Region II consists of modes that are outside the current hori- 
zon, but were inside the horizon at the time of the interaction, r v . Using the relevant 
expansions for the mode functions pi) we find the momentum space propagator 



F{l{k^T 2 ) 



XH 2 (A; 2 r 1 r 2 ) 



(27T) 



XH 2 (k 2 nT 2 y 



(27T) 



Ak 3 



Ak A 

hj 



sm(2kr v — ne) 



T, 



v 

00 ds 



l-(AlRT v ) 2 * 1 //X 
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sin(2s — ne) 



m 

2e 
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sin(2s 



where in the second line we changed variables s = —kr v , and as the main contribution to the 
integral comes from the lower limit, we moved the upper limit to infinity. Written thus, the 
two remaining integrals match to the leading order in e, yielding / = ^§ sin(2s — ire) w 
0.063 + 0(e). Converting back to position space we then have 



F 1L (x 1 ,x 2 




l+2e 



A2e 



2e 

2(A /fi r ) 



2c 



log 



r 2 
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This result should be contrasted with the infrared contribution (14). Both asymptote to 
a constant value for very late times, but the IR contribution is dominant due to an extra 
factor of e _1 . One can again verify that the apparent divergence as r 2 — > with t% fixed is 
due to the assumption T\ ~ r 2 , and vanishes as one computes the potentially divergent part 
more carefully. 

The UV modes: Finally, in region III the scale k is inside the current horizon between 
t and r 2 . Using the expansion ^ the contribution to the propagator ^ becomes 



*lt(*,Ti,7s) 



XH 2 Tir 2 f T2 dr v . 

— cos k{Ti - t v ) sin fc(r 2 - t v ) 



(2tt) 2 2k 2 
XH 2 nr 2 
(2tt) 2 2k 
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(15) 



where we again defined s = — kr v and moved the upper limit to infinity. To leading order in 
e we then find 



XH 2 T X T 2 



1 - (A/ R r 2 ) 2e 
2e 



2 \H 



(2tt) 2 2A; 

• (cos fc(n + r 2 )Ci(-2A;r 2 ) - sinA;(ri + r 2 ) (| + Si(2A;r 2 ; 
where the Cosine and Sine integrals are defined as 



CifaO 



cos s 



ds. 



and Si(x) 



Transforming back to position space we then get 
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ds. 
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2e 2 



rf/c fe (cos fc(n + r 2 )Ci(2A;r 2 ) - sin/c(ri + r 2 ) (| + Si(2A;T 2 )j) 
- (A/iiT 2 ) 2e 



+ 



£V 



2e 2 

Above, we noted that the upper limit of the integral doesn't contribute due to the oscillatory 
nature of the integrand, and then obtained a rough estimate of the magnitude of the integral 
by the value of the integrand at the lower limit multiplied by frequency of oscillation. We 
again find a contribution that asymptotes to a constant for very late times, but is sub- 
dominant to the infrared and intermediate contributions found above. 
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F. The one-loop propagators 



We can now write down the propagators to one loop using the results above. The tree 
level contributions from the infrared modes ([7]) to the propagators can be written as 



F QL (x 1 ,x 2 ) 



rl3U ZJ2 /•min(-r 1 1 ~r 2 1 ) 

e *(^) F(A;jTl ^ )ra _ (Tl75|) e^ dkk ~^ (16) 
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(17) 



Adding the dominant one loop IR contribution (14) to the tree contribution Fol(^i, ^2)1 we 
find the full propagator for late times as 



F(xx,x 2 ) 



H 2 1 
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n <"> T"2, 



(18) 



where the symmetrization T\ <->■ r 2 corresponds to including the mirror images of the dia- 
grams in figure [3j This shows that for a massive field (where e > 0) the comoving infrared 
cutoff Am drops out of the late time correlators, which match those obtained from a Bunch- 
Davies/Euclidean vacuum computation. 



G. Early time propagator 

We should contrast the propagator found above with the corresponding result for early 
times or very light fields. The infrared contribution to the tree level propagator for light (or 



massless) fields can be computed from (16) by a power series expansion in e as 



F^(x 1 ,x 2 ) « ^ (log(A IR r ) + log ( min( ^' r2) )) (1 + 0{e)) , (19) 

where the first logarithm measures how far the infrared cutoff is from the initial horizon, 
and the second logarithm gives the number of e-folds of inflation. The one-loop correction 
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can be similarly computed by expanding ( 13 ) as a power series in e. After some algebra this 
gives 



F( L {x 1 ,x 2 ) 



XH 2 
9(2tt) 4 
1, r 2 
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ii 
8 Vif 
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31og(A //? r ) 



12 \H 



0{e)) . 



(20) 



These results showcase the familiar logarithmic divergences for late times r 2 — >■ 0, but as 
shown earlier, these are absent for truly late times if e > 0. The solution to this apparent 



paradox is that the expansions in (19), (20) are valid when e C 1/ log {tq/t^) , which means 
the logarithmic expansion breaks down for sufficiently late times, and the propagator is 



better described by (18), which is well behaved. 



III. DISCUSSION AND CONCLUSIONS 



Massive scalar We conclude from section II F that in the distant future, the loop cor- 
rections are insensitive to the choice of infrared cutoff. The result for the comoving cutoff 
approaches that of the de Sitter invariant Euclidean vacuum results as r| e , i.e. as a negative 
power of the scale factor. 



Massless scalar, or a light scalar at early times Here the result of section |II G| is ap- 
plicable. This can be the relevant behavior during the entire course of slow roll inflation. 
The difference between the result of section |II G| and section II F provides a measure of 



the theoretical uncertainty in the predictions of the two-point function caused by uncer- 
tainties in the initial state at the start of inflation. This difference is proportion to a 
A log 2 (T2/V0) log (AjflTo) relative to the tree level contribution. The slow roll conditions for 
pure A0 4 give the constraint that the initial expectation value for must be larger than the 
Planck scale. The condition that the density fluctuations today be sufficiently small (10~ 5 ) 
gives the constraint that A ~ 10~ 15 . So we see in this example the effect will be at least 
10 -11 smaller than the tree-level piece, even after 60 e-folds of inflation. We conclude there- 
fore, that while these secular logarithms are present during primordial slow-roll inflation, 
matching the magnitude of our observed density perturbations constrains the parameters of 
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the inflaton model such that the coefficient of the secular log is very small, unless inflation 
is very long-lived. 

Massless scalar at late times The Euclidean vacuum methods break down in this case 
(see for example [33j for a recent discussion of this situation) due to the nonexistence of a de 
Sitter invariant vacuum for a massless scalar |34| . However the comoving IR cutoff approach 
still provides a well-defined perturbative expansion, up until the point when A log 3 (T2/70) ~ 



1 (see (20)). The prediction is that the theoretical uncertainty in correlators becomes large 
at sufficiently late times. This is the practical sense in which de Sitter spacetime suffers a 
quantum instability. In the slow roll example, this does not happen until about 10 5 e-folds, 
so is really only of academic interest as far as primordial inflation goes. 

Graviton at late times However it should be noted our result for the massless scalar 
matches qualitatively with the computation of |23j . There the full one- loop correction for 
graviton self-energy was computed with the same kind of comoving infrared cutoff that we 
have used in the present paper. Compared to the tree-level result the contribution is of 
order (^jf~^ l°g 2 ( r 2/ r o)- Again this is very small over the 60 e-folds needed for primordial 
inflation. However over the present epoch such a term may become important, if we are 
indeed evolving toward a period where dark energy dominates. To apply these ideas to our 
present epoch, we could imagine simply shifting tq to the present time, and placing a comov- 
ing infrared cutoff just outside our present horizon. This suggests quantum uncertainties 
due to the indeterminancy of our initial state would only become large after (^-jj^^j ~ 10 122 
e-folds of late time expansion, where now H w 10~ 33 eV is the Hubble parameter today. 
It is interesting that this provides a mechanism for the possible breakdown of semiclassical 
physics in a de Sitter region, as predicted by quantum gravity arguments in [30J. 

Acknowledgments 

This research is supported in part by DOE grant DE-FG02-91ER40688-Task A and an 
FxQI grant. We thank Richard Easther, Richard Woodard and Helmut Friedrich for helpful 



15 



discussions. 



[1] H. Friedrich, "On The Existence Of N-Geodesically Complete Or Future Complete Solutions 

Of Einstein Field-Equations With Smooth Asymptotic Structure", Communications In 

Mathematical Physics 107 (1986), no. 4, 587-609. 
[2] H. Friedrich, "Existence and structure of past asymptotically simple solutions of einstein's 

field equations with positive cosmological constant", Journal of Geometry and Physics 3 

(1986), no. 1, 101 - 117. 
[3] H. Friedrich, "On the global existence and the asymptotic behavior of solutions to the 

Einstein-Maxwell- Yang-Mills equations", J.Diff.Geom. 34 (1991) 275-345. 
[4] J. H. Traschen and C. Hill, "Instability Of De Sitter Space On Short Time Scales", Phys.Rev. 

D33 (1986) 3519-3525. 

[5] N. Tsamis and R. Woodard, "Relaxing the cosmological constant", Phys.Lett. B301 (1993) 
351-357. 

[6] R. H. Brandenberger, H. Feldman, and V. F. Mukhanov, "Classical and quantum theory of 
perturbations in inflationary universe models", |arXiv:astro-ph/9307016. 



[7] N. Tsamis and R. Woodard, "Strong infrared effects in quantum gravity", Annals Phys. 238 
(1995) 1-82. 

[8] N. Tsamis and R. Woodard, "The Quantum gravitational back reaction on inflation", Annals 

Phys. 253 (1997) 1-54, |arXiv:hep-ph/96023"l6l 
[9] N. Tsamis and R. Woodard, "Quantum gravity slows inflation", Nucl.Phys. B474 (1996) 



235-248, arXiv:hep-ph/9602315 



[10] V. F. Mukhanov, L. W. Abramo, and R. H. Brandenberger, "On the Back reaction problem 
for gravitational perturbations", Phys. Rev. Lett. 78 (1997) 1624-1627, arXiv:gr-qc7960 9026 

[11] L. W. Abramo, R. H. Brandenberger, and V. F. Mukhanov, "The Energy - momentum tensor 
for cosmological perturbations", Phys.Rev. D56 (1997) 3248-3257, [arXiv:gr-qc/9704037| 



[12] L. Abramo and R. Woodard, "One loop back reaction on chaotic inflation", Phys.Rev. D60 



(1999) 044010, |aIXiv:astro-ph /9811430 



[13] S. Weinberg, "Quantum contributions to cosmological correlations", Phys.Rev. D72 (2005) 



043514, arXiv:hep-th/0506236 



16 



[14] S. Weinberg, "Quantum contributions to cosmological correlations. II. Can these corrections 



[15 
[16 

[17 



[18 

[19 

[20 

[21 
[22 

[23 

[24 

[25 

[26 

[27 
[28 
[29 



become large?", Phys.Rev. D74 (2006) 023508, arXiv:hep-th/0605244 



A. Polyakov, "Beyond space-time", |arXiv:hep-tn /0602011 

M. van der Meulen and J. Smit, "Classical approximation to quantum cosmological 



correlations", JCAP 0711 (2007) 023, arXiv:0707.0842 



N. Bartolo, S. Matarrese, M. Pietroni, A. Riotto, and D. Seery, "On the Physical Significance 
of Infra-red Corrections to Inflationary Observables", JCAP 0801 (2008) 015, 
larXiv:0711.4263l 

A. Polyakov, "De Sitter space and eternity", Nucl.Phys. B797 (2008) 199-217, 
larXiv:0709.2"899l 

P. Adshead, R. Easther, and E. A. Lim, "The 'in-in' Formalism and Cosmological 



Perturbations", Phys.Rev. D80 (2009) 083521, arXiv:0904.4207 



C. Burgess, L. Leblond, R. Holman, and S. Shandera, "Super-Hubble de Sitter Fluctuations 



and the Dynamical RG", JCAP 1003 (2010) 033, arXiv:0912.1608 



A. Polyakov, "Decay of Vacuum Energy", Nucl.Phys. B834 (2010) 316-329, arXiv:0912.5503 



D. Krotov and A. M. Polyakov, "Infrared Sensitivity of Unstable Vacua", Nucl.Phys. B849 



(2011) 410-432, arXiv:1012.2107 



N. Tsamis and R. Woodard, "One loop graviton selfenergy in a locally de Sitter background", 



Phys.Rev. D54 (1996) 2621-2639, arXiv:hep-ph/9602317 



D. Seery, "Infrared effects in inflationary correlation functions", Class. Quant. Grav. 27 (2010) 
124005, |arXiv: 1005. 1649 



E. Kahya, V. Onemli, and R. Woodard, "A Completely Regular Quantum Stress Tensor with 
w < -1", Phys.Rev. D81 (2010) 023508, |a7Xiv:0904.4811 



S. Kundu, "Inflation with General Initial Conditions for Scalar Perturbations", 
larXiv:1110.4688l 

E. W. Kolb and M. S. Turner, "The Early universe", Front.Phys. 69 (1990) 1-547. 



D. H. Lyth, "The curvature perturbation in a box", JCAP 0712 (2007) 016, arXiv:0707.0361 



D. A. Lowe, "Some comments on embedding inflation in the AdS/CFT correspondence", 



Phys.Rev. D77 (2008) 066003, arXiv:0710.3564 



[30] D. A. Lowe and S. Roy, "Punctuated eternal inflation via AdS/CFT", Phys.Rev. D82 (2010) 



063508, arXiv: 1004. 1402 



17 



[31] K. Enqvist, S. Nurmi, D. Podolsky, and G. Rigopoulos, "On the divergences of inflationary 
superhorizon perturbations", JCAP 0804 (2008) 025, |arXiv:0802. 0395 



[32 
[33 
[34 
[35 

[36 

[37 

[38 

[39 

[40 

[41 

[42 

[43 

[44 
[45 

[46 

[47 

[48 



J. Maciejko, "An introduction to non-equilibrium many-body theory", 2007, unpublished. 



S. Hollands, "Massless interacting quantum fields in deSitter spacetime", arXiv: 1105. 1996 



B. Allen, "Vacuum States in de Sitter Space", Phys. Rev. D32 (1985) 3136. 

S. L. Adler, "Massless, Euclidean quantum electrodynamics on the five- dimensional unit 

hypersphere", Phys. Rev. D6 (1972) 3445-3461. 

S. L. Adler, "Massless Electrodynamics On The Five-Dimensional Unit Hypersphere: An 
Amplitude - Integral Formulation", Phys. Rev. D8 (1973) 2400-2418. 
I. T. Drummond, "Dimensional Regularization of Massless Theories in Spherical 
Space-Time", Nucl. Phys. B94 (1975) 115. 

I. T. Drummond and G. M. Shore, "Conformal Anomalies for Interacting Scalar Fields in 
Curved Space-Time", Phys. Rev. D19 (1979) 1134. 

I. T. Drummond and G. M. Shore, "Dimensional Regularization of Massless Quantum 
Electrodynamics in Spherical Space- Time. 1", Ann. Phys. 117 (1979) 89. 

D. Marolf and I. A. Morrison, "The IR stability of de Sitter QFT: results at all orders", Phys. 



Rev. D84 (2011) 044040, arXiv:1010.5327 



D. Marolf and I. A. Morrison, "The IR stability of de Sitter: Loop corrections to scalar 



propagators", Phys. Rev. D82 (2010) 105032, arXiv: 1006.0035 



S. Hollands, "Correlators, Feynman diagrams, and quantum no-hair in deSitter spacetime", 
larXiv: 1010.53671 

D. Marolf and I. A. Morrison, "The IR stability of de Sitter QFT: Physical initial conditions", 
larXiv:1104.4343l 

A. A. Starobinsky, "Stochastic De Sitter (Inflationary) Stage In The Early Universe", 1986. 
A. A. Starobinsky and J. Yokoyama, "Equilibrium state of a selfinteracting scalar field in the 



De Sitter background", Phys. Rev. D50 (1994) 6357-6368, |arXiv:astro -ph/9407016 

S. R. Coleman and F. De Luccia, "Gravitational Effects on and of Vacuum Decay", Phys. Rev. 

D21 (1980) 3305. 

K. Goldstein and D. A. Lowe, "Initial state effects on the cosmic microwave background and 



trans-planckian physics", Phys. Rev. D67 (2003) 063502, hep-th/ 0208167. 



K. Goldstein and D. A. Lowe, "Real-time perturbation theory in de Sitter space", Phys. Rev. 

18 



D69 (2004) 023507, hep-th/0308135 

[49] K. Goldstein and D. A. Lowe, "A note on alpha-vacua and interacting field theory in de 
Sitter space", Nucl. Phys. B669 (2003) 325-340, hep-th/0302050l 

[50] For conformally coupled matter, this was studied long ago |35| 136] . and for more general 
matter in |37H39| . These questions have been revisited more recently in |33[ |4"DT|43| . 

[51] The landscape of string theory provides a wealth of viable vacua typically separated by 
Planck scale potential walls. This fits well with the picture of eternal inflation, where a 
Planck scale curvature foam tends to dominate the volume on any given space like slice. It is 
worth emphasizing that approximate approaches to performing computations in these 
situations, such as stochastic inflation |44[ 145] . where super horizon fluctuations in a fixed de 
Sitter background are considered in order to compute an initial wavefunction for the 
universe, or Coleman-de Luccia bubble nucleation [46J, where tunneling to an open universe 
sets the initial conditions for inflation, both require extensive fine tuning. In |44t 145] . because 
only the dynamics of superhorizon modes is considered, the high curvature spacetime foam is 
ignored, and replaced by a fixed de Sitter background, which is an inconsistent 
approximation. Likewise if we wish to apply the results of [46j to eternal inflation in the 
string landscape we must assume we are in a special state where the semiclassical instanton 
methods are applicable. The philosophy in the present paper is to sidestep these issues, by 
making minimal assumptions about the initial state of inflation, in order to better model this 
generic strongly curved foam outside the initial homogeneous region. 

[52] It should be noted that a more general family of de Sitter invariant vacua |34j . known as the 
alpha vacua, are possible for quantum fields on a fixed de Sitter background. See for example 
|47H49] and references therein. However this approach is likely to break down when quantum 
gravity fluctuations are included. This is essentially because the Euclidean/Bunch-Davies 
vacuum is the only one that preserves the desired short distance singularities of correlators, 
whereas the alpha vacuum correlators contain additional singularities at antipodal points, 
creating light-like non-locality. This in turn can lead to strong back-reaction effects once 
spacetime fluctuations are included. 



19 



